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Underground muon events detected by the MACRO experiment at Gran Sasso have been studied for different 
purposes. The studies include the vertical muon intensity measurement, multiplicity distribution, lateral and 
angular muon distribution and searches for substructures inside muon bundles. These analyses have contributed 
to bring new insights in cosmic ray physics, in particular in the framework of primary cosmic ray composition 
studies. Moreover, this activity allows the testing and tuning of Monte Carlo simulations, in particular for aspects 
associated with models of hadronic interactions and muon propagation through the rock. 



1. INTRODUCTION 

Muons detected deep underground are a useful 
tool for different physics and astrophysics items. 
These muons are the decay product of mesons 
originating in the very first hadronic interac- 
tions in the atmosphere or in secondary interac- 
tions during the shower development. Therefore, 
their study can provide many informations about 
primary cosmic rays (CR) composition and/or 
high energy hadronic interactions. Muons ar- 
riving in the underground Gran Sasso Labora- 
tory have crossed at least h > 3100 hg/crri^ of 
standard rock, corresponding to a muon energy 
cut Ef^ > 1.3 TeV at the surface. This means 
that the primary CR energies range from some 
TeV/nucleon up to the maximum energies well 
beyond the "knee" . 

The MACRO experiment has collected a large 
amount of muon data in the last decade, at a rate 
of ~ 6.6 X 10® muon events/live year, and the ~ 
6% of these are multiple muon events. Many un- 
derground observables have been studied. The 
multiplicity distribution, i.e. the rate of muon 
events as a function of their multiplicity, is a 
quantity strongly dependent on the primary com- 
position model. A detailed analysis on primary 
composition has been performed by the MACRO 
collaboration [0,^: one of the result is that data 
prefers a composition model with an average mass 



slightly increasing with energy above the knee[|. 

Nevertheless, in the context of composition 
studies the knowledge of the hadronic interaction 
model is crucial. The main contribution to the 
systematic uncertainties in these analyses are due 
to the interaction models adopted and it is impor- 
tant to find out new observables to test the relia- 
bility of these models implemented in the Monte 
Carlo simulations. Moreover, is intrinsically im- 
portant to test interaction models in kinematical 
regions not yet explored at accelerators or collid- 
ers. For instance, about 5% of MACRO muon 
data comes from pp interactions with ^/s > 2 
TeV, while about 30% of the muons observed are 
the decay products of mesons produced with a 
(pseudo)rapidity r]cm > 5. The situation is more 
evident if we consider that part of primary in- 
teractions in the atmosphere are nucleus-nucleus 
interactions, and very few data for energy Eiab ^ 
150 A GeV are available. In the following, we will 
focus on the decoherence analysis and on cluster 
analysis, two different tools able to extract infor- 
mations on the interaction model adopted in the 
simulation codes. 

2. THE MACRO DETECTOR 

The MACRO detector Q is a large area appa- 
ratus located in hall B of the Gran Sasso Lab- 
oratory at an average depth of 3800 hg/cm^ of 
standard rock (E'^ > 1.3 TeV). It has a mod- 
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ular structure, organized in six almost identi- 
cal "supermodules" covering an horizontal sur- 
face area of ~ 1000 m^. The apparatus is 
equipped with three different and independent 
sub-detectors: streamer tube chambers for par- 
ticle tracking, scintillator counters for timing and 
energy loss reconstruction and nuclear track-etch 
detectors optimized for the search of magnetic 
monopoles. The wire view of the streamer tube 
system is complemented with a second view, dis- 
posed at 26.5° with respect to the wire view, real- 
ized with aluminium pick-up strips. The spatial 
resolutions are crn/=l.l cm and 175'= 1.6 cm for the 
wire and strip view respectively. This arrange- 
ment allows the 3-D track reconstruction with an 
intrinsic pointing resolution of 0.2°. 

3. SIMULATION TOOLS 

The full simulation chain used to interpret 
MACRO data is composed of an event genera- 
tor modelling high energy hadronic interactions, 
included in a shower propagation code which fol- 
lows particles above threshold up to a given at- 
mospheric depth. A muon transport code prop- 
agate muons inside the mountain overburden the 
apparatus and a detector simulator produces an 
output in the same format of real data. 

• Event generator 

Most of the analyses of MACRO data have 
been carried on using the original HEMAS inter- 
action model based on the parametrizations 
of minimum bias events collected at the SppS col- 
lider at CERN Q . According to these results, the 
charged multiplicity is sampled from a negative 
binomial distribution and the transverse momen- 
tum contains a power law component. The model 
includes nuclear target effects and extrapolations 
to higher energies are performed in the context of 
log(s) physics. 

The interaction model in Ref. Q (called 
"NIM85" ) has been tested in a parametrized form 
with MACRO muon data |8|jl|l. This model ne- 
glects some important experimental results: for 
instance, the charged multiplicity is sampled from 
a Poisson distribution and the transverse momen- 
tum distribution contains only a pure exponential 



functional form. 

Presently new interaction models are under 
study: DPMJET |l§, QGSJET [|ll, SIBYLL 
|l^ and HDPM, the original interaction model 
of the shower simulation code CORSIKA [T^ . 
These are phenomenological models where the in- 
teractions are treated at parton level, with the ex- 
ception of the HDPM generator, which is based 
on the DPM model but it is built according to 
parametrized results. They have the common 
feature to refer to the Regge-Gribov theories for 
the modelling of the soft part of the interactions, 
where perturbative QCD cannot be applied. Nev- 
ertheless, the transverse component of the inter- 
actions is not constrained by the theory and is 
introduced "by hand" , according to some experi- 
mental results such as the seagull effect or the 
Cronin effect |Q in nuclear interactions. In this 
context, it is useful the comparison of the model 
one to each other to estimate the systematic un- 
certainty associated to the unknown transverse 
structure of the interactions. 

• Cascade code 

Two different shower propagation codes have 
been used in the analyses: HEMAS0 and 
CORSIKA |l|. HEMAS is conceived as a fast 
tool to compute the hadronic, electromagnetic 
and muon component of air showers. It has 
been extensively used in the MACRO collabora- 
tion and has revised many improvements since its 
first release. It introduces some approximations 
in the shower development, so that the model can 
be used to follow only particles with a minimum 
energy in atmosphere E > 500 GeV. The elec- 
tromagnetic size is computed by means a semi- 
analytical method. It can be used interfaced with 
the DPMJET and HEMAS interaction models. 

Recently, the CORSIKA @ Monte Carlo 
code, generally used in surface EAS-arrays, has 
been interfaced with the muon transport code to 
propagate muons up to the Gran Sasso depth. 
This code has been used only in the analysis of 
high multiplicity events. 



^the name HEMAS refers both to the hadronic interaction 
model and to the cascade code 



• Muon transport code 

The muon propagation in the rock has been 
reahzed using the PROPMU package , which 
represents an improvement with respect to the 
propagation model included in the original 
HEMAS version. It takes into account muon 
energy loss due to multiple Coulomb scattering 
and to discrete processes, such as bremsstrahlung, 
pair production processes and photonuclcar inter- 
actions. 

• Detector simulator 

The response of the apparatus is simulated 
by means a GEANT based code, called 
GMACRO. It reproduces all relevant physical 
processes in the detector and produces an output 
in the same format of real data, so as to process 
data using the same offline chain of real data. 

4. ANALYSIS RESULTS 

4.1. Decoherence Function 

The decoherence function, defined as the dis- 
tribution of the muon pair separation in multiple 
muon events, is mainly connected with the muon 
lateral distribution with respect to the shower 
axis and is very sensitive to the transverse struc- 
ture of the hadronic interaction models. The de- 
coherence distribution is instead weakly depen- 
dent on the primary composition model, as far as 
the shape of the distribution is concerned. There- 
fore, the study of this function allows to some 
extent the disentangle between the two effects. 
MACRO has studied the decoherence function up 
to the maximum distance allowed by the appara- 
tus (~ 70 m). The unfolded decoherence func- 
tion (the distribution corrected for the detector 
effects) is shown in FigJ^, where it is compared 
with the predictions of the HEMAS Monte Carlo. 

A further check on the reliability of the sim- 
ulation code has been performed in different 
rock/zenith windows, constraining some compo- 
nent of the shower development with respect to 
others Again, the comparison shows a good 
agreement between data and Monte Carlo. 

Finally, the study of the decoherence function 
in the very low distance region has shown that 
the QED process fi^ + N ^ + N + fi+ + fi^ 
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Figure 1. Experimental decoherence distribu- 
tion (corrected for the detector effects) compared 
with the predictions of the HEMAS Monte Carlo 
(MACRO-fit primary composition model) 



at small distances must be taken into account if 
we want to reproduce the experimental data [pT| . 
This is shown in Fig.^, where the low distance re- 
gion of the decoherence function is shown before 
and after the correction. The contribution of this 
process is negligible compared to the e+e~ pair 
production process in the GeV range, but it be- 
comes progressively more important in the TeV 
region p^ ]. 

4.2. Cluster Analysis 

The search for substructures inside muon bun- 
dles is able to provide additional information with 
respect to traditional methods jl^. In some 
events muon bundles appear to be splitted into 
"clusters" and we ask if this feature is the re- 
sult of simple statistical fluctuations in the muon 
lateral distribution or if there is some dynami- 
cal correlation connected with the development 
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(HEMAS and NIM85). Most of the effect can be 
explained as the result of fluctuations of the muon 
density inside the bundles. Considering that the 
density (namely the average number of muons 
per unity of area) depends both on primary mass 
number and on the modelling of meson transverse 
momentum, we can expect that the cluster ef- 
fect is sensitive to the composition model and to 
the hadronic interaction model at the same time. 
On the other hand, if we consider the interaction 
models quoted in previous sections, the sensitiv- 
ity of the cluster effect on the interaction model 
becomes weaker. This is shown in Fig.^ where is 
reported the cluster rates as a function of the pa- 
rameter Xcut for different interaction models and 
for a fixed primary composition model (the model 
derived from the fit of the MACRO multiplicity 



Figure 2. Zoom of the low distance region of the 
decoherence function before and after the inclu- 
sion of muon pair production process, and com- 
parison with the Monte Carlo. In the inset it is 
shown the difference of the distributions before 
and after the correction, compared with the ex- 
pectations of the Monte Carlo. 



of the shower in the atmosphere. From an ex- 
perimental point of view, we select muon bundles 
with at least 8 muons underground {N^ > 8) cor- 
responding to CR primary energies Eprimary ~ 
1000 TeV. The search for muon clusters is per- 
formed by means of an iterative cluster finding 
algorithm, which groups the muons depending on 
the choice of a free parameter called Xcut (for the 
definition of this parameter see Ref. |p^ ). 

In Ref. ||l^ has been pointed out that this 
method is sensitive both on the primary com- 
position model and on the hadronic composi- 
tion model. The comparison was made between 
two extreme composition model (the "heavy" and 
"light" composition models ||2^) and between 
two very different hadronic interaction model 



N-Cluster Event Fractions in a projected view 
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Figure 3. Fraction of events with 1,2 and 3 clus- 
ters in muon bundles for different hadronic in- 
teraction models. The detector effects have been 
considered and all the simulations have been ob- 
tained with the MACRO-fit composition model. 
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distribution |2|). In this case we are forced to en- 
hance the sensitivity applying some selection cri- 
teria which correlate the underground substruc- 
tures with the hadronic interaction features in the 
atmosphere. 

Apart from statistical fluctuations, a Monte 
Carlo study has revealed other two mechanisms 
responsible for the cluster effect: 

• muons belonging to the same cluster have a 
larger probability to have a common parent me- 
son in the steps of shower tree generation; 

• muons belonging to the same cluster are the 
decay products of mesons highly correlated in the 
phase space of the very first hadronic interactions 
in the atmosphere. We considered only events 
reconstructed as two-cluster events by the algo- 
rithm with a fixed Xcut and we studied the kine- 
matical variables pseudorapidity 77 and azimuthal 
angle (j) of the parent mesons in the first interac- 
tion of the shower. The typical topology of these 
selected events is a muon rich cluster close to the 
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shower axis, generally the remnant of the shower 
development after many steps, and a far cluster 
with few muons directly generated in the very 
first hadronic interactions. Fig.^ shows the dis- 
tributions of the relative distance in the rj — cj) 
space of pairs of muon parent mesons originated 
in the first interactions: the topological selection 
of muons belonging to the same or to different 
clusters reflects in a strong selection of different 
phase space regions of hadronic interactions. 

A quantitative computation of the relative con- 
tributions of these effects is at presently under 
study. 

5. CONCLUSIONS 

The MACRO experiment has collected a large 
amount of muon data in the last decade. The 
analyses performed with these data allowed to 
draw conclusions about several items of muon 
physics. The dimension and granularity of the 
detector allowed the detection of multiplicity up 
to ^ 40, under more than 3000 hg/cm^ of rock 
overburden. The study of the multiplicity distri- 
bution has been used to extract informations on 
the primary cosmic ray composition. 

The modelling of the transverse component of 
the interaction models at TeV energies, connected 
with the lateral distribution of cosmic ray shower, 
is one of the main source of uncertainties in the 
simulation codes. The analysis of the decoherence 
function has shown the reliability of the HEMAS 
Monte Carlo, used in most of the MACRO anal- 
yses. 

The study of second order effects, like the 
search for jet substructures inside muon bundles, 
is a useful tool to add new informations with re- 
spect to conventional analyses. The physical in- 
terpretation of the results has shown the dynam- 
ical origin of the effect, connected with the devel- 
opment of the shower in the atmosphere. 



Figure 4. 2-dimensional distributions of the dis- 
tance of first interaction mesons in the rj—c/) space. 
It is shown the case of no selection, of the selection 
of muons belonging to the same cluster (poorest 
or richest) or to different clusters. 
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